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Abstract. We discuss the evolution of quantum correlations for a system of two 
two-level atoms interacting with a common reservoir. The Markovian master equation 
I is used to describe the evolution of various measures of quantum correlations. It is 

shown that different measures of quantum correlations exhibit qualitatively different 
behaviour in their evolution. 
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1. Introduction 

It is now well established that quantum entanglement, which is a measure of 
quantum correlations, is a necessary resource to perform some quantum information 
algorithms [Ij. However, it has been realized that quantum entanglement does not 
include all quantum correlations, and there are separable states which exhibit quantum 
correlations different from entanglement. It became important to recognise the nature of 
correlations and properly discriminate quantum correlations from the classical ones. A 
number of measures of quantum correlations has been introduced, and the most popular 
of them is quantum discord [2l|3] (see also the review [1] and papers cited there). It is 
crucial, from the point of view of future application, to know how quantum correlations 
evolve in time when a multipartite system interacts with the dissipative environment. 

The time evolution of entanglement for a system of two qubits, or two two-level 
atoms, can be qualitatively described for various physical situations, and it has been 
studied extensively in recent years [514T2]. A lot of discussion has been devoted to the 
problem of disentanglement of the two-qubit system in a finite time, despite the fact 
that all the matrix elements of the two-atom system decay only asymptotically. Yu and 
Eberly [9] coined the name "entanglement sudden death" (ESD) to the process of finite- 
time disentanglement. ESD has recently been confirmed experimentally |Tp]. Another 
problem related to the entanglement evolution that attracted attention is the evolution 
of the entangled qubits interacting with the non-Markovian reservoirs [T314T5] . It has 
also been shown that squeezed reservoir leads to the steady-state entanglement [IE] and 
revivals of entanglement [T7] . 

Beside entanglement, a considerable amount of attention has been paid in recent 
years to quantum discord. Generally, quantum discord is difficult to calculate because 
it requires the optimisation procedure which usually amounts to extensive numerical 
calculations. In case of two qubits analytical results can be obtained for some 
specific families of states, like Bell-diagonal states or states having maximally-mixed 
marginals [18]. Another interesting family of states is a family of the two-qubit X-states, 
which are of interest here. Ali et al. [19] reported a closed form solution for quantum 
discord of the X states. However, it turned out that their algorithm is not universal. 
Lu et al. [20] have proven that the universal set of orthogonal projective measurements 
cannot be found for the full family of X states. Some counterexamples have been given 
in [20l[2Tj . Chen et al. [21], however, confirmed applicability of the algorithm for several 
special cases of X states. Lu et al. have found that the probability distribution of 
measurements is centralised around a specific von Neumann measurement, which they 
called the maximal-correlation-direction measurement. This observation justifies the 
Ali et al. algorithm in a statistical sense, because, it turns out that for 99.4% of cases of 
numerically generated X states the algorithm gives correct results. The situation with 
arbitrary two-qubit states is more complicated and the best what has been achieved so 
far are the two transcendental equations obtained by Girolami and Adesso [22j which 
must be solved numerically to get quantum discord. 
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Dakic et al. ^B] have introduced another measure of quantum correlation based on 
the Hilbert-Schmidt distance measure, which is called geometric discord. The advantage 
of this measure is that it allows for getting analytical formulas for general two-qubit 
states. Girolami and Adesso [2l] introduced the so called observable measure of bipartite 
quantum correlation which is a lower bound to the geometric discord. Behaviour of the 
geometric discord under decoherence has been studied in [25]. Recently Bellomo et 
al. [26] compared the dynamics of geometric and entropic quantifiers of the different 
kinds of correlations in a non-Markovian open two-qubit system under local dephasing. 

In this paper we study the evolution of concurrence (a measure of entanglement), 
quantum discord and geometric discord in a system of two two-level atoms interacting 
with a common reservoir being in a vacuum state. The evolution of the system 
is described by the Markovian master equation introduced by Lehmberg [27] and 
Agarwal [28], taking into account the cooperative behaviour of the atoms. It is shown 
that different measures of quantum correlations evolve in time quite differently. We 
compare their evolution for a family of pure initial states. Since the evolution of the 
system is described by a realistic master equation, which is a good testing ground for 
studying physical processes involving two atoms, we believe the results obtained shed 
new light on understanding quantum correlations. 

2. Master equation 

We consider a system of two two-level atoms, A and i?, with ground states \gi) and 
excited states |ej) {i = A,B) connected by dipole transition moments /Xj. The atoms 
are located at fixed positions and and coupled to all modes of the electromagnetic 
field being in the vacuum state. 

The reduced two-atom density matrix evolves in time according to the Markovian 
master equation given by [271 - 129] 



where (5") are the raising (lowering) operators, and S- is the energy operator of 
the ith atom, Tu = T are the spontaneous decay rates. We assume that the two atoms 
are identical. The parameters Tij and flij {i 7^ j) depend on the distance between the 
atoms and describe the collective damping and the dipole-dipole interaction defined, 
respectively, by 



dp 
di 



2 2 



i^uji[S^,p]-iY^ Qij [S^Sj , p] 



i=l i^j 




(1) 
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and 

^ 3^ / cos kra sin /en,- cos/cnA , . 

. = -r ^ + 1 + ^ , (3) 

where k = ojq/c, and r^j is the distance between the atoms. Here, we assume, with no 
loss of generahty, that the atomic dipole moments are parallel to each other and are 
polarised in the direction perpendicular to the interatomic axis. 

To describe the evolution of the two-qubit system the standard basis of atomic 
product states can be used: |1) = le^) ® le^), |2) = |e^) ® Igs), |3) = Iqa) ® \^b), 
|4) = \gA)®\gB)- It is easier, however, to find the solutions of the master equations when 
using instead of the standard basis, a basis of the collective states: |e) = |e^)®|eB), \s) = 

75 (Ica) ® Iqb) + \9b) ® les)), \a) = 75 (le^) ® Igs) - \9a) ® Icb)), \g) = Iqa) ® Iqb)- 
The states |s) and \a) are the symmetric and antisymmetric states of the two-atom 
system. They are maximally entangled states, or Bell states of the two-atom system. 
Assuming that initially the system density matrix has the so called X form, which is 
preserved during the evolution according to the master equation ([1]), and we get the 
following system of equations for the density matrix elements [2S] 

Peeit) = pee(0)e-2r*, 

Pss{t) = p..(0)e-^+* + Pee(O)^ [e-^+* - e-^^*] , 

Paait) = PaaiO)e-^-' + Pee(O)^ [c"^"* - C'^^*] , (4) 

Pasit) = p..(0)e-(™-)*, 
Peait) = p,e(0)e-(r+2^-»)*, 

where F-t = F ± Tab- For calculating the quantum correlations measures we need the 
solutions for the density matrix elements in the standard product basis, which can be 
expressed in terms of the matrix elements ([1]) in the following way 

Pll{t) = Pee{t), 
Pliif) = Pegif), 

P22{t) = [Pss{t) + Paa{t) + Pas{t) + Psa{t)]/2, (5) 
P33(t) = [Pssit) + Paa(t) + Pas(t) - Psa(t)]/2, 
P23it) = [pssit) - paa{t) + Pas{t) " psa{t)]/2. 

The solutions (jl]) and ^ are used to find the evolution of various measures of quantum 
correlations. 

3. Measures of quantum correlations 

3.1. Entanglement 

Probably the most celebrated and studied manifestation of quantum correlations is 
quantum entanglement. To quantify the entanglement, various entanglement measures 
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Figure 1. Evolution of concurrence (a), quantum discord (b), and geometric quantum 
discord (c), for the two-photon Bell-like state (fT3| and the interatomic distance 
rAB = A/8. The white lines show section of the surface for p = 2/3. 

have been introduced. We use here concurrence introduced by Wootters [30]. In the 
case of X states we consider, the concurrence can be calculated analytically, and it has 
the form [7] 

C{t) = max{0,Ci(t),C2(t)} 

Ci(t) = 2 (IpmWI - Vp22{t)p33{t)) (6) 

C2{t) = 2 (\p2s{t)\ - ^pu{t)p,,{t)) 

Inserting into ([6]) the solutions (jl]) and (j5]), we find the values of Ciit) and C2(t), 
and whenever one of the two quantities becomes positive, there is a some degree of 
entanglement in the system. 
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Figure 2. Evolution of correlations for the state p3p and interatomic distance 
TAB = A/8 for p = 2/3. The evolution is split into two time regions (notice the 
different scales of both parts). 



3.2. Quantum discord 

To calculate quantum discord we use the algorithm introduced by Ali et al. [19], which 
works quite well in our case. We have checked for many numerically generated X states, 
for which the algorithm fails to find correct extremum, that the differences between the 
numerically found extrema and the values given by the algorithm are so small that they 
cannot be resolved in the scale of the figures. In our calculation, we assume that the 
measurement is performed on the subsystem B. For X state we have 

V =5(p^)-5(p) + min{/Ci,/C2}, 



^2= - ^ pa logs Pa - S{p^) 



(7) 



where >S'(p), S{p^) mean the von Neumann entropy for the two- atom system and the 
subsystem B, respectively, sf = Tr[p(T^] = pu + p22 — P33 — P44 is the third component 
of the Bloch vector of the subsystem A, and rj = {pul + \p23\- H{x) is the Shannon 



entropy. 



Quantum correlations in a two- atom system 



7 



0.2 



0.15 



0.1 



0.05 




N ^ N 



1 

Tt 



0.05 r 







C 


0.04 




V 






g 






G 


0.03 








Figure 3. Same as Fig. [31 but for p = 1, i.e. initially both atoms excited. 



3.3. Geometric quantum discord 

Geometric quantum discord is defined as 

g = 2mm\\p- (8) 

where A* is a set of zero-discord states, and 1 1 . . . 1 1 denotes tlie Hilbert-Sclimidt norm. 
Tfie factor 2 in front has been introduced for normahsation. The general formula for 
the geometric quantum discord takes the form 

?B| |2 _|_ 1 12^1 12 



.) (9) 

where s^ is the Bloch vector for the subsystem B (we assume that measurement is 
performed on subsystem _B), and k^ax is the largest eigenvalue of the matrix 

ir = s^(s^)* + T*T, (10) 

where the superscript t means transposition. 

For X states we have 

Q = 

Ql = 4(|pi4| 

Q2 = 2(|/0i4p - IP2I 

where sf = pn — P22 + P33 — P44 is the third component of the Bloch vector for the 
subsystem B and = TT[p{af <S> cx^), {i,j = x,y,z), are elements of the correlation 
matrix. 



min{^i,^2} 

|2 



IP23 



2\2 



(11) 
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The observable measure of Girolami and Adesso [2^ is given by the formula 



G = ^ 2Tr(K) - ^J6tT{K^) -2[Tt{K)]^ 



(12) 



with the matrix K given by ( !T0|) . 

4. Evolution of quantum correlations 

To compare the evolution of various measures of quantum correlations in a system of 
two two-level atoms governed by the master equation ([1]), we assume that the initial 
state is the Bell-like state 



where p is the population of the upper state |1). In Fig. [T]we illustrate the evolution 
of concurrence, quantum discord and geometric quantum discord for the whole range 
of values of p. For p < 1 the state (fT3|) is a superposition of the states |1) and |4) 
which has nonzero two- photon coherence pi4, while for p = 1 it is a product state of 
both atoms excited. It is seen from Fig. [1] that the evolution of different measures of 
quantum correlations depends on the value of p and shows essential differences, not 
only quantitative but also qualitative. The white line seen on the figures indicates the 
section of the surface at p = 2/3. The evolution for this value of p is illustrated in more 
details in Fig. O It has been shown before [H] that there is sudden death and revival 
of entanglement in the system for p > 0.5, which is clearly seen in Fig. [2)d. The other 
correlations do not exhibit sudden death, but there is an increase in quantum discord 
for rt > 2. The most striking feature is the behaviour of geometric quantum discord 
with a deep crater for Ft < 1 as seen in Fig. [T]3, and the cusp seen in Fig. [2^. The 
observable measure of geometric quantum discord given by (fT2l) and shown in Fig. [2^ 
smoothes out the cusp of geometric discord. 

For p = 1 the state f|T3l) becomes the product state of both atoms being excited for 
which all the correlations are zero. However, as it is evident from Fig. [H concurrence 
remains zero for some time, while the other correlations increase immediately after 
the start of the evolution. It is clearly evident from Fig. [3l For short times the 
quantum discord increases, reaches the maximum and goes down to the minimum in 
order to increase again to the subsequent maximum, and eventually it goes down to zero 
asymptotically. This behaviour has been reported in [32]. The concurrence remains zero 
up to Ft ~ 4 and next abruptly becomes nonzero, the effect that has been referred to as 
sudden birth of entanglement [12]. Both revival of entanglement and the sudden birth 
of entanglement are due to collective behaviour of the two atoms when the interatomic 
distance is smaller than the wavelength of light emitted by individual atom [TT | [T2 l [33] . 
Here we take the interatomic distance equal to A/8. 

From Fig. [3^ it is seen that the geometric quantum discord has a minimum with 
sharp cusps, and the observable measure of geometric quantum discord, which is a lower 
bound for geometric quantum discord, is a really tight bound, except for the interval 



^) = v^ii) + v^r^|4) 



(13) 
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around the minimum, where the bound is not so tight. For this time interval geometric 
quantum discord behaves quite differently than quantum discord, and both differ from 
concurrence. 

5. Conclusions 

We have discussed the dynamics of various measures of quantum correlations in a two- 
atom system interacting with a common reservoir in a vacuum state. The evolution 
of the system is described by the Lehmberg-Agarwal Markovian master equation, 
which takes into account collective behaviour of the atoms. The collective spontaneous 
emission is a source of quantum correlations in the system. We have shown that the 
evolution of different measures of quantum correlations is qualitatively different, with 
a rather strange behaviour of the geometric discord. Some aspects of the evolution of 
quantum correlation in such a system has been studied in [Mj. Recently Piani [35] 
has argued that the geometric discord is not a good measure for the quantumness of 
correlations. 
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